非耦合、弱耦合正倒向随机微分方程的高阶数值方法及误差估计

【摘要】1990年,Pardoux和Peng首次证明了非线性倒向随机微分方程(简称,BSDE)解的存在性和唯一性[46].在[49],Peng首先得到了BSDE和抛物型偏微分方程(PDE)的关系,且在[48]中研究了基于BSDE的最优控制问题的随机最大值原理,随之BSDE的相关理论和数值计算分析发展的十分迅速.相应地,国际上正倒向随机微分方程(简称,FBSDEs)也随着正向随机微分方程(简称,SDE)和BSDE的发展得到了广泛深入的理论和应用研究.而且FBSDEs的研究成果在许多领域,如金融数学、风险度量、随机控制、博弈理论、物理、化学、航天以及随机偏微分方程等数学领域有很重要的应用[4,9-11,17,19-21,23,26,27,32,37,38,40,42,44,45].至今,FBSDEs理论已经成为概率论及随机分析理论中的一个重要分支.但是到目前为止,FBSDEs的计算理论还并不丰富,绝大多数计算结果的精度较低,很多这方面问题的研究尚处于初级阶段,难以达到实际应用的要求.随着FBSDEs科学计算的深入研究,需要研究求解FBSDEs的新方法.结合FBSDEs理论和科学计算理论,需要用创新性思维提出求解FBSDEs的新方法,研究方法的收敛性、稳定性,进行理论数值分析和误差估计.这些研究对深刻理解FBSDEs的特性、丰富发展FBSDEs科学计算理论、加快FBSDEs实际应用都有重要的意义。本文旨在从FBSDEs的理论出发,研究FBSDEs的高精度数值求解方法,理论误差估计以及算法实现.到目前为止,关于数值求解FBSDEs的数值方法大概分为两类：一类是基于FBSDEs和相应抛物型PDE的关系,通过求解PDE的粘性解达到求解FBSDEs的目的；另一类是通过对FBSDEs直接数值离散得到相应的数值计算格式来求解FBSDEs的数值解.我们研究的求解FBSDEs的高精度数值格式的方法是后一类方法.本论文主要从正倒向随机微分方程的特性出发,结合确定性科学计算理论,系统研究非耦合正倒向随机微分方程(简称,decoupledFBSDEs)、弱耦合正倒向随机微分方程(简称,weaklycoupledFBSDEs)以及带跳的非耦合正倒向随机微分方程(简称,decoupledFBSDEswithjumps)的高阶数值求解方法,严格地数值理论分析所提格式的稳定性、收敛性,得到格式的误差估计,通过数值模拟分析格式的高效高精度性,研究所提方法在金融、随机控制等方面的应用。下面将介绍论文的主要内容以及论文的结构.第一章引言,介绍了第三章到第六章中我们研究的主要问题.第二章,介绍了正向随机微分方程和带跳的正向随机微分方程解的存在唯一性定理,列举一些简单的求解正向随机微分方程和带跳的正向随机微分方程的数值格式以及真解和数值解之间的关系.第三章,主要研究非耦合正倒向随机微分方程(decoupledFBSDEs)的一步数值计算方法.我们研究的方法是：(1)利用随机分析和FBSDEs理论,通过条件数学期望,将FBSDEs变为相对确定的积分方程；(2)利用科学计算中的数值积分理论,在时间方向上用数值积分高精度近似得到的确定性方程中的积分；(3)用Gauss数值积分公式准确近似(2)中的条件期望；(4)用数值代数,求解(3)得到的方程组,得到高精度格式；(5)利用随机分析理论、FBSDEs理论和科学计算理论,分析格式的稳定性、收敛性和误差.本章的创新之处：通过随机分析、FBSDEs理论和确定型科学计算理论的结合,提出求解FBSDEs的高精度数值格式.具体创新点如下：(1)构造一种新的特殊的标准布朗运动；(2)利用构造的新的布朗运动,提出求解非耦合正倒向随机微分方程的高精度数值格式,得到稳定性、收敛性和误差估计；(3)对求解非耦合正倒向随机微分方程的高精度数值格式进行算法实现.本章来自于：Zhao,W.,Zhang,W.,andJu,L.,ANumericalMethodanditsErrorEstimatesfortheDecoupledForward-BackwardStochasticDifferentialEquations,Commun.Comput.Phys.15(2014),pp.618-646.在第三章研究内容的启示下,第四章主要研究求解非耦合正倒向随机微分方程的多步数值计算方法.首先我们提出全新的倒向正交多项式的概念,基于该倒向正交多项式构造新的特殊的标准布朗运动；然后利用构造的新的布朗运动,提出非耦合正倒向随机微分方程的高精度数值格式,严格地理论证明了这一新的多步数值格式的稳定性,得到了多步方法的高阶误差估计；最后通过数值实验同样得到了高精度数值格式的收敛阶,从而验证了理论结果的正确性.本章来自于：Zhao,W.,Zhang,W.,andJu,L.,Amultistepschemefordecoupledforward-backwardstochasticdifferentialequations,已投稿.第五章,考虑弱耦合正倒向随机微分方程的相关问题,其中弱耦合是指正向随机微分方程的系数函数是依赖倒向随机微分方程的解Y的,但不依赖于过程Z.对该问题进行研究得到了高精度的数值格式.(1)改进了一种求解弱耦合正倒向随机微分方程的数值格式,这种格式我们称之为改进的Euler数值格式和改进的Markovian迭代数值格式.随后对原数值格式和改进的数值格式进行了严格的理论误差分析,时间方向上证明了原Euler数值格式和改进的Euler数值格式都具有一阶的时间收敛性,原Markovian迭代数值格式和改进的Markovian迭代数值格式均是收敛的.(2)利用第三章中使用的方法对(1)所提出的改进的数值格式进行了时间方向上的收敛性证明,证明了所提格式具有一阶的时间离散高精度收敛性,数值算例验证了这一占创新之处：(1)改进的Euler数值格式和改进的Markovian迭代数值格式在计算时明显减少了求解条件期望需要的计算量,在时间上得到了极大的改善.(2)收敛性分析表明,收敛精度从原来的半阶提高到一阶,同样是另一个极大的改盖本章来自于Zhang,W.,Zhao,W.,Euler-typeschemesforweaklycoupledforward-backwards-tochasticdifferentialequationsandtheoptimalconvergenceanalysis,FrontiersofMath-ematicsinChina,DOI:10.1007/s11464-014-0366-6,2014.Errorestimatesofthenumericalschemeforsolvingweaklycoupledforward-backwardstochasticdifferentialequations,已完成.第六章,研究带跳的正倒向随机微分方程的相关理论,提出求解带跳的正倒向随机微分方程的高精度数值格式,并进行相关的理论误差分析.首先,明确了补偿泊松过程驱动的带跳的正倒向随机微分方程和补偿泊松随机测度驱动的带跳的正倒向随机微分方程的不同以及各自的性质.然后,考虑由补偿泊松随机测度驱动的带跳的倒向随机微分方程和带跳的非耦合正倒向随机微分方程的问题.研究带跳的正向随机微分方程的伊藤泰勒展开公式,随后提出新的算子应用于求解的方程,给出求解问题方程的参考方程.接着根据参考方程提出求解问题方程的高精度数值格式,并对所提的数值格式进行严格的理论误差分析,最后用数值算例来验证所提格式的高阶收敛性.本章来自于：Anewkindofaccuratediscrete-timeapproximationofbakcwardstochasticdiifferen-tialequationwithjumps.已完成.Secondordernumericalschemeanderrorestimatesofsolvingdecoupledforwardbackwardstochasticdifferentialequationswithjumps.己完成.以下是本文的章节目录及主要结论.一、第一章引言；二、第二章正向SDE和带跳的正向SDE；三、第三章非耦合FBSDEs的一步数值格式和误差估计；四、第四章非耦合FBSDEs的多步数值格式和误差估计；五、第五章弱耦合FBSDEs的改进数值格式及其误差估计；六、第六章带跳的非耦合FBSDEs的高精度数值格式及误差估计.第三章：首先构造一种新的标准的布朗运动△Ws,利用这种新的布朗运动提出一种新的非耦合正倒向随机微分方程的一步高精度数值格式.理论上,在一般情形下严格地分析了所提格式的误差,然后当正向随机微分方程使用不同精度的格式求解时,对每一种情况给出了特殊的误差分析；特殊地,当正向随机微分方程使用弱2.0阶数值格式进行求解时,所提出的求解正倒向随机微分方程的数值格式是二阶高精度格式.随后给出数值算例,表明所提格式的精度的确是二阶的,同时也验证了理论误差结论是正确的.考虑下面的非耦合正倒向随机微分方程：通过引入新的布朗运动：得到求解decoupledFBSDEs的一步数值格式：格式3.1.1.(Xn,Yn,Zn)是(Xt,Yt,Zt)在时间t=tn处的数值解,通过下面的方程组来求解(Xn,Yn,Zn),其中n=N-1,N-2,…,1,0,这里φ表示的是经典的求解SDE的数值格式,例如Euler格式、Milstein格式以及一些强收敛格式、弱收敛格式等等.对上面的一步数值格式,我们有下面的误差结论.定理3.3.1.令(Xt,Yt,Zt),t∈[0,T]以及(Xn,Yn,Zn),n=0,1,…,N,分别表示非耦合FBSDEs(3.1)的精确解和由格式3.1.1得到的数值解.假设函数f(t,X,Y,Z)对X,Y和Z满足Lipschitz连续性,并且Lipschitz常数是L.这里c0是由(3.2)定义的时间剖分的正则性系数.那么对于足够小的时间步长Δtn,下面的不等式成立这里n=N-1,….,1,0,其中C是一个依赖于c0和L的正常数,C’是一个依赖于co,T和L的正常数Riy和Riz在方程(3.9)和(3.16)中有定义,并且引理3.4.1.如果f(t,x,y,z)∈C2,4,4,4b,b(t,x),σ(t,x)∈C2,4b,φ∈C4+α,α=∈(0,1)和|b(t,x)|2≤K(1+|x|2),|σ(t:x)|2≤K(1+|x|2),那么对充分小的时间步长Δt。,我们对任意0≤n≤N-1,有其中C是一个依赖于T,K以及系数函数b,σ,f和终端函数φ导数上确界的正常数.引理3.4.2.引理3.4.1的条件成立,那么对充分小的时间步长Δt。,我们对任意0≤n≤N-1,有其中C是一个依赖于T,K,以及系数函数b,σ,f和终端函数φ导数上确界的正常数.引理3.4.3.引理3.4.1的条件成立,那么对充分小的时间步长Δt。,我们对任意0≤n≤N-1,有其中C是一个依赖于T,K,以及系数函数b,σ,f和终端函数φ导数上确界的正常数.引理3.4.4.引理3.4.1的条件成立,那么对于充分小的时间步长Atn,对任意的0≤n≤N-1,有其中C是一个依赖于T,K,以及系数函数b,σ,f和终端函数φ导数上确界的正常数.定理3.4.1.假设2.3.2和引理3.4.1的条件成立.那么对于充分小的时间步长At。,我们对任意0≤n≤N-1,有其中C1是一个依赖于c0,T和L的正常数,C2是一个依赖于C0,T,L,K,方程(3.1)的解Xt的初值,以及系数函数b,σ,f和终端函数φ导数上确界的正常数.第四章：首先从概念上提出了一类新的正交多项式{Pi(s)}Li=0,我们称之为倒向正交多项式.随后研究了提出的倒向正交多项式的性质.利用该倒向正交多项式,构造了一种新的布朗运动△Wtn,8=∫stnPKz(r)dWr,结合这一新的布朗运动和倒向随机微分方程,提出了求解正倒向随机微分方程的新的多步数值格式.严格地理论证明了这一新的多步数值格式的稳定性,得到了多步数值格式的高阶误差估计.最后通过求解数值算例,验证了理论误差分析和数值计算结果一致.在这一章,仍然考虑第三章的非耦合正倒向随机微分方程,在本部分将给出多步数值格式,定义和主要的结论.定义4.1.1.我们称一组定义在区间[0,1]上的多项式{Qi(s)}Li=0为倒向正交多项式,如果Qi(s)∈Pi[0,1],i=0,1,…,L,并且下面的条件成立引理4.1.1.定义在区间[0,1]上的多项式集合{Qi(s)}Li=0满足条件(4.3).令Qi(s)=引理4.1.2.令{Pi(s)}Li=0是区间[a,b]上的倒向正交多项式.那么我们有P0(α)=1和Pi(α)>1,1≤i≤L.引理4.1.3.令{Pi(s)}Li=0是区间[a,b]上的倒向正交多项式.假设函数g(s)的直到i+1阶导数g(i+1)有界,那么对0≤i≤L,我们有格式4.2.1.令K=max[Ky,Kz,Kf).{Xn)Nn=0是正向随机微分方程的数值解.给定随机变量X0,YN-i和ZN-i,i=0.1,…,K-1那么对n=N-K,…,0,1.使用下面的格式解Zn：2.使用下面的格式解Yn：引理4.3.2.假设N和K是两个非负整数并且N≥K,△t是任意正数.假设{an},{bn)和{Rn},n=0,1,…,N是非负的,并且满足下面的不等式其中C是一个依赖于C1,C2,C3,K和T的正常数.定理4.3.1.令(Xt,Yt,Zt),t∈[0,T]和(Xn,Yn,Zn),n=0,1…,N分别表示非耦合FBSDEs的真解和由多步数值格式4.2.1得到的近似解.假设函数f(t,Xt,Yt,Zt)对Xt,Yt和Zt是Lipschitz连续的,并且Lipschitz常数是L.令c0是由(4.6)定义的时间剖分的正则参数.令K=max{Ky,K2,Kf),B=max{Bky,BKf}以及PKz=其中C是一个依赖于c0,T,L,B,K,QKz(0)和PKz的常数Riy和Riz分别由(4.14)和(4.19)定义,同时引理4.4.2.令Rng和Rnz是定义在参考方程(4.14)和(4.19)中的局部截断误差.那么在假设4.4.1下,我们有下面的局部估计：其中C≥0是一个依赖于T,以及系数函数b,σ,f和终端函数φ导数上确界的常数.定理4.4.1.假设4.4.1和假设2.3.2成立,并且假设初始值满足Mey=O((△t)Ky+1)其中C是一个依赖于c0,T,L,K,B},QKz(0),PKz,方程(4.1)的解Xt的初始值,以及系数函数b,σ,f和终端函数φ导数上确界的常数.第五章：在这一章节,介绍了一种新的求解弱耦合正倒向随机微分方程的Euler数值格式和它的迭代格式—Markovian迭代数值格式.尽管改进的数值格式和原来出现在[6]中的数值格式有一些共性,但是改进的格式计算量明显比原来的数值格式小.随后对改进的Euelr数值格式进行了误差分析和理论证明,证明了改进的Euler数值格式具有一阶收敛性,从而改进了原来的求解弱耦合正倒向随机微分方程的Euler数值格式的半阶收敛性结论,最后给出了数值算例同样验证了格式的收敛精度是一阶的.在这一部分我们考虑弱耦合正倒向随机微分方程：格式5.1.3.(改进的Euler数值格式)给出随机变量X0,YN和ZN.对n=N-1,N-2,…,0,我们使用下面的数值格式来求解随机变量Xn,Yn和Zn格式5.1.4.(改进Euler数值格式的Markovian迭代格式)给出X0,m,YN,m=φ(XN，m)并且令un,0(X0,m)=0.对n=N-1,N-2,…,0,我们使用下面的数值格式来求解Xn,m,Yn,m和Zn,m定理5.3.1.(i)令(Xt,Yt,Zt),t∈[0,T]和(Xn,Yn:Zn),n=0,1,…,N,分别是弱耦合FBSDEs方程(5.1)的精确解和格式5.1.1或者格式5.1.3的数值解.假设5.2.3和假设5.2.4成立.那么对于充分小的时间步长△t,对n=N-1,N-2,…,0,下面的估计式成立(ii)进一步,如果(Xn,m,Yn,m,Zn,m),n=0,1,…,N和m=0,1,2,…,是格式5.1.2或者格式5.1.4的近似解,那么存在常数C>0以及0<c<1使得,对充分小的时间步长Δt,我们有下面的估计式成立这里C是一个不依赖于时间剖分和数值格式近似解的正常数.引理5.3.2.在假设5.2.4下,对充分小的时间步长Δt,我们对任意0≤n≤N-1有下面的误差估计：其中C是一个不依赖于时间剖分和近似解的正常数.引理5.3.4.令(Xt,Yt,Zt),t∈[0,T]以及(Xn,Yn,Zn),n=0,1,…,N,是弱耦合FBSDEs(5.1)的精确解和由方程(5.17),(5.18)和(5.19)获得的近似解.假设函数f(t,X,Y,Z)对X,Y和Z是一致Lipschitz连续的.那么对于充分小的时间步长At,对n=N-1,N-2,…,0,下面的估计式成立其中Riy和Riz是由(5.14)和(5.15)定义的,并且并且如果假设5.2.4成立,我们有下面的估计：使用(5.24),有这里C是一个不依赖与时间剖分和方程(5.17)-(5.19)解的正常数.定理5.3.2.令(Xt,Yt,Zt),t∈[0,T]和(Xn,Yn,Zn),n=0,1,…N分别是弱耦合FBSDEs(5.1)的精确解和格式5.1.3的近似解.假设5.2.3和假设5.2.4成立.那么对充分小的时间步长At,对任意0≤n≤N-1,下面的估计式成立|Ytn(x)-Yn(x)|2≤C(1+|x|2)(△t)2,其中C是一个正常数.第六章：本章研究了带跳的正向随机微分方程的伊藤泰勒扩展公式,随后提出新的标准布朗运动和新的泊松跳运动,利用这两种新的运动给出问题的参考方程,用于提出求解带跳的倒向随机微分方程和带跳的非耦合正倒向随机微分方程的高精度数值格式,并对所提的数值格式进行严格地理论误差分析,最后用数值计算来验证所提格式的高阶收敛性.本章主要研究如下由补偿泊松随机测度驱动的带跳的非耦合正倒向随机微分方程：其中ξ=φ(XT),Γ:=∫Eρ(e)U(e)λ(de),sup|ρ(e)|≤K.这里W是一个d-维布朗运动,μ是一个独立于w的补偿泊松测度,并且β(de,ds)=μ(de,ds)-λ(de)ds.在第三章构造的新的布朗运动的启发下,为了高精度的求解r,大胆的猜测可以构造类似的补偿泊松测度基于此,给出求解(X,Y,Z,Γ)的高精度数值格式：格式6.3.1.给出随机变量X0,YN,ZN和TN.在定义(6.11)和定义(6.17)中,令s=(?)n+1我们得到△Wtn+1和△(?)tn+1(0≤n≤N-1)的定义.对n=N-1,N-2,…,0,通过下面的方程组求解随机变量Yn,Zn和Γn并且Xn+1=Xn+Φ(tn,Xn,△tn,△Wn,1,ξn+1).例如：Xn+1可以由下面的弱二阶格式解出引理6.4.2.如果f(t,x,,y,z)∈C2,4,4,4b,b(t,x),σ(t,x)∈C2,4b,φ∈C4+α,n∈(0,1)并且假设2.4.1成立,那么对充分小的时间步长Δtn,我们对任意0≤n≤N-1,有E[|Rny|2]≤C(1+E[|Xn|8])(△t)6,其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.3.如果引理6.4.2的条件成立,那么对充分小的时间步长Δtn,任意0≤n≤N-1,有其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.4.如果引理6.4.2的条件成立,那么对充分小的时间步长Δtn,任意0≤n≤N-1,有其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.5.如果引理6.4.2的条件成立,那么对充分小的时间步长△tn,任意0≤n≤N-1,有|Rnz|2≤C(1+EXntn[|xn|8])(△t)6,其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.6.如果引理6.4.2的条件成立,那么对充分小的时间步长Δtn,任意0≤n≤N-1,有其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.7.如果引理6.4.2的条件成立,那么对充分小的时间步长Δtn,任意0≤n≤N-1,有其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.引理6.4.8.如果引理6.4.2的条件成立,那么对充分小的时间步长Δtn,任意0≤n≤N-1,有|RnΓ|2≤C(1+EXntn[|Xn|8])(△t)6,其中C是一个依赖于T,K和系数函数b,σ,f以及终端函数φ导数上界的正常数.定理6.4.1.令(Xt,Yt,Zt,Γt),t∈[0,T]和(Xn,Yn,Zn,Γn),n=0,1,…,N,分别是非耦合正倒向随机微分方程(6.1)的精确解和由数值格式6.3.1求得的数值解.假设函数f(t,Xt,Y,Zt,Γt)对Xt,Yt,Zt和гt是Lipschitz连续的,并且Lipschitz常数是L.令c0是由定义(6.4)给出的时间剖分正则参数.那么对充分小的时间步长Δt,下面的估计式成立其中n=N-1,…,1,0,C是一个依赖于c0和L的常数,C’是一个依赖于c0,T,L.由(6.9)和(6.16)定义的截断误差Riy和Riz,以及求解正向随机微分方程带来的弱收敛误差项Rny1,Rny2,Rnz1,Rnz2,RnΓ1和RnΓ2的常数,他们分别定义为：定理6.4.2.假设2.4.2和引理6.4.2的条件成立.那么对于充分小的时间步长△tn,任意0≤n≤N-1.有其中C是依赖于c0和L的常数,C1是一个依赖于c0,T和L的常数,C2是依赖于c0,T,L,K,方程(6.1)中Xt初始值,以及系数函数b,σ,f和终端函数φ导数上界的常数.
【作者】张微；
【导师】赵卫东；
【作者基本信息】山东大学，概率论与数理统计，2014，博士
【关键词】非耦合正倒向随机微分方程；弱耦合正倒向随机微分方程；带跳的倒向随机微分方程；带跳的非耦合正倒向随机微分方程高精度；数值计算格式；误差估计；

【参考文献】
[1]康彬.高温PEMFC不锈钢双极板的电化学性能研究[D].大连海事大学,材料工程,2014,硕士.
[2]路文焕.高校数字化教育资源的法律保护[D].天津大学,教育经济与管理,2004,硕士.
[3]何小娟.前列地尔对体外培养类缺血再灌注模型视网膜神经节细胞的影响[D].广西医科大学,眼科学,2013,硕士.
[4]刘坤.河北省蔬菜价格波动特征及影响因素研究[D].河北经贸大学,产业经济学,2013,硕士.
[5]马冉冉.中等职业学校学生英语学习动机研究[D].首都师范大学,学科教学,2013,硕士.
[6]潘勇.吉林省农村公路养护内容及政策研究[D].长安大学,交通运输工程,2013,硕士.
[7]曾尤,王函.纳米炭/环氧树脂复合材料黏弹阻尼性能[J].科学通报,2014,33:3232-3239.
[8]魏敬亮.信息化逐级而上[J].中国石油企业,2004,11:80-81.
[9]付长民,底青云,王妙月.海洋可控源电磁法三维数值模拟[J].石油地球物理勘探,2009,03:358-363+386+255.
[10]李鹏飞.HPV38 E6/E7对角质形成细胞自噬的影响及机制探讨[D].第四军医大学,皮肤病与性病学,2014,硕士.
[11]郑鸿飞.CO_2在油饱和多孔介质中的扩散实验研究[D].大连理工大学,能源与环境工程,2013,硕士.
[12]欧娟.主流叙事与另类话语[D].河南大学,中国现当代文学,2004,硕士.
[13]李娜,马静,孙书利.带多丢包和滞后随机不确定系统的最优线性估计[J].自动化学报.
[14]梁怀阳.十二醇/三聚氰胺—甲醛/脲醛树脂微胶囊相变材料的制备与表征[D].厦门大学,化学工程,2014,硕士.
[15]肖海涛.气调包装生鲜牛肉腐败菌生长模型的构建[D].上海海洋大学,食品工程,2013,硕士.
[16]刘志强.给水输配系统优化设计的研究[D].天津大学,环境工程,2004,硕士.
[17]杨志勇.流动人口基本公共服务均等化研究[D].内蒙古大学,公共管理,2014,硕士.
[18]李青.莱辛《三四五区间的联姻》的主题研究[D].西安外国语大学,比较文学与世界文学,2013,硕士.
[19]吴建德.基于神经网络的负载模拟器控制方法研究[D].西北工业大学,机械电子工程,2004,硕士.
[20]李攀峰.清华大学继续教育的品牌定位与传播研究[D].电子科技大学,工商管理（专业学位）,2012,硕士.
[21]邱青松.转FSHα/β基因克隆猪分子鉴定及常规生理生化指标测定[D].江西农业大学,动物遗传育种与繁殖,2012,硕士.
[22]刘清泉.我国石油重力勘探的回顾和展望[J].石油地球物理勘探,1988,06:719-727+765.
[23]桂婷婷.论我国公司机会规则的完善[D].复旦大学,民商法学,2012,硕士.
[24]焦大蒙.加工中心链式刀库及机械手可靠性试验方法研究[D].吉林大学,机械工程,2014,硕士.
[25]周鹏.汪元量《醉歌》、《湖州歌》、《越州歌》叙事研究[D].宁夏大学,中国古代文学,2014,硕士.
[26]刘秀秀.不同钾效率类型油菜干物质积累及钾素吸收利用差异研究[D].华中农业大学,植物营养学,2014,硕士.
[27]刘铮峰.广西钦州坭兴陶溯源及发展的初步研究[D].广西民族大学,科学技术史,2013,硕士.
[28]吕伟.烟草企业FMS中激光导引式AGV防撞保护功能与优化路径调度的设计与实现[D].浙江工业大学,2012.
[29]王蕾.杜氏藻中β-胡萝卜素及甾醇的研究[D].内蒙古大学,化学,2014,硕士.
[30]邵良杉,付曙光,薛立军.企业核心竞争力的Web挖掘研究[J].计算机系统应用,2005,08:91-94.
[31]聂纯阳.基于碳纳米管及其复合物电极的电容去离子技术研究[D].华东师范大学,2013.
[32]李桂兴.力学知识的综合应用[J].中学生理科月刊,2004,05:44-48.
[33]张春生.PCGC编译码算法的研究及FPGA实现[D].中南民族大学,通信与信息系统,2013,硕士.
[34]王帅.影视资源在中学语文课程教学中的开发与应用策略研究[D].东北师范大学,广播电视艺术学,2012,硕士.
[35]林丹.CRM对企业运作模式的影响[D].厦门大学,2001.
[36]陈康.永磁无刷直流电动机控制方案的研究[D].南京航空航天大学,电机与电器,2013,硕士.
[37]罗薇.全面成本管理在建筑工程中的应用[D].西南交通大学,结构工程,2013,硕士.
[38]谢玲燕.小学语文阅读教学中关注言语形式的实践与研究[D].华中师范大学,教育管理,2012,硕士.
[39]李新江.面向视频监控的自动人数统计研究[D].西南交通大学,信号与信息处理,2014,硕士.
[40]孙婉华.高空抛物致人损害案件法律问题研究[D].沈阳师范大学,法律,2013,硕士.
[41]王銮.随机时滞系统的控制问题研究[D].鲁东大学,运筹学与控制论,2014,硕士.
[42]万福磊.建筑结构连续性倒塌数值模拟方法研究[D].中国建筑科学研究院,2011.
[43]蔡俊青.造纸专用碳酸钙制备及其动力学研究[D].河北工业大学,化学工艺,2004,硕士.
[44]马超.激光熔覆镍基复合涂层组织与性能研究[D].大连理工大学,材料无损检测与评价,2013,硕士.
[45]任翠霞.TC4表面激光诱导原位制备TiC/TiB_x钛基复合涂层显微组织及性能的研究[D].东北大学,材料学,2011,硕士.
[46]翟蕾.Web文档分类及信息抽取的研究与开发[D].西北工业大学,2001.
[47]钱驰.基于静电纺丝的可控喷印技术和组织工程应用[D].苏州大学,控制理论与控制工程,2013,硕士.
[48]张亚飞.中学生信息共享行为研究[D].河北大学,情报学,2014,硕士.
[49]漆永新.解读钢铁企业信息化[J].中国钢铁业,2004,02:5-8.
[50]黄毅,陈永胜.石墨烯的功能化及其相关应用[J].中国科学(B辑:化学),2009,09:887-896.